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1 Introduction 



An interesting question is how present and future experiments will be able to probe the couplings of the 
Higgs boson at a high level of precision, see Ref. [1] for a discussion. There is a wide variety of beyond 
the Standard Model (BSM) theories where the Higgs couplings differs from the Standard Model (SM) ones 
by less that 10%, as discussed in Ref. [2]. Among many papers dealing with the subject we quote those 
in Refs. [3-7]. For the most recent update on the subject we refer to the work of Refs. [8-10]. Interim 
recommendations to explore the coupling structure of a Higgs-like particle can be found in Ref. [11]. 

In this work, following Independence Day [12,13], we imagine that there is a huge space of theories, 
represented by local and renormalizable Lagrangians where SM is only one point. A possible strategy to 
look for deviations from the SM is the following: 

• we take the SM as the theory of "light" degrees of freedom, i.e. d = 4 operators 

• we simulate the unknown extension of the SM by the most general set of d = 6 operators that are 
obtained by integrating out the heavy degrees of freedom (we also assume no sensitivity to operators 
with d > 8 at LHC). This is equivalent to say that the BSM theory is unknown or matching is too 
difficult to carry out, so we write the most general set of interactions consistent with symmetries. The 
effective theory contains an infinite number of operators but only a finite number is needed for present 
(LHC) precision. 

With enough statistics it should be possible to fit a,-, the Wilson coefficients of the d = 6 operators, and there 
are two possibilities: a) they are close to zero (where zero = SM) or b) they are not. Option a) tell us that 
NLO corrections (or the residual theoretical uncertainty at NNLO level) and the a,- coefficients are small, 
and the SM is actually a minimum in our Lagrangian space or very close to it. This will explain nothing but 
it is internally consistent. Option b) raises serious problems since the effect of local operators is large and 
they cannot be included only at LO, but inserting operators in SM loops creates even more problems. 

In case it is option b) we should move in the Lagrangian space and adopt a new renormalizable La- 
grangian with the virtue of making zero that specific (large) Wilson coefficient a, ; local operators are then 
redefined w.r.t. the new Lagrangian. Of course there will be more Lagrangians projecting into the same set 
of operators but still we could see how our new choice handles the rest of the data. 

In principle, there will be a blurred arrow in our space of Lagrangians, and we should simply focus 
the arrow. Without invoking the explicit example of Supersymmetry this is the so-called inverse problem 
introduced in Refs. [14,15]: if LHC finds evidence for physics beyond the SM, how can one determine the 
underlying theory? 

It is worth noting that this question is highly difficult to receive a complete answer at the LHC. The main 
goal will be to identify the structure of the effective Lagrangian (i.e. the different scalings of the various 
d = 6 operators) and to derive qualitative information on new physics; the question of the ultraviolet (UV) 
completion cannot be answered unless there is sensitivity to d > 6 operators. Therefore, we are looking for 
a relatively modest goal on the road to understand if the effective theory can be UV completed (bottom-up 
approach with no obvious embedding). 

To set up our definitions of an enlarged theory we have to specify the concept of Higgs fields: it is the set 
of scalar fields that break electroweak symmetry (EW) by developing a vacuum expectation value (VEV). 
What we are looking for is evidence of SM Higgs properties or deviations from the SM behavior; in the 
latter case one has to understand consistency with other EW symmetry-breaking frameworks. Alternatively 
we can consider scenarios with more scalar fields, that are not Higgs fields (Higgs partners); the problem 
with more VEVs, or one VEV different from (T , Y) = (± , 1) (T is isospin and Y is hypercharge), is partially 
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related to the rho-parameter [16] which at tree-level is given by 



1 Li Ci | v t | 2 +r,« 2 



= 7}(7} + l)-if, 



r i = T i {T i +\) 



Plo = x 



2 L^ 2 |v ; - 



(1) 



where the sum is over all Higgs fields, v ; -(w,) gives the VEV of a complex(real) Higgs field with hypercharge 
Yi and weak-isospin 7}. Our considerations will be presented in Section 6.5. The experimental limit on 
p — 1 are rather stringent. For a complete discussion of models respecting custodial symmetry we refer to 
Ref. [17]. 

In this paper we do not discuss questions related to spin, mass or CP quantum numbers but only cou- 
plings. In particular we discuss couplings to vector bosons since they control the unitarity behavior of 
longitudinal VV -scattering at high energy [18] (automatic in the SM). We also discuss the effects of Higgs 
partners and of Higgs self-couplings in Section 8. 

For a better illustration of our approach we observe that a consistent effective theory, defined by 



has arbitrary Wilson coefficients a" which, however, give the leading amplitudes in an exactly unitary S- 
matrix at energies far below A. The theory is non-renormalizable, which means that an infinite number of 
higher operators must be included. Nevertheless there is a consistent expansion of amplitudes in power of 
E/A. Our goal will be to understand the d = 6 operators as a first step towards an UV completion, possibly 
a weakly-coupled one i.e. one where weakly-coupled new physics opens up around A and restores unitarity 
(a different scenario, classicalization, has been proposed in Refs. [19,20]). 

In other words, the question is: can we classify the low-energy (LHC) observables that determine the 
road to UV completion? Note that there is a claim in the literature [21,22] that the coefficients a,- must be 
positive to have an UV completion which respects the usual axioms of S -matrix theory. In particular, in 
the work of Ref. [22] it is shown that UV completion is encoded in the sign of the scattering amplitude for 
longitudinal vector-bosons and that weakly-coupled UV completion requires a positive sign. Constraints 
on the sign of the couplings in an effective Higgs Lagrangian using prime principles have been derived in 



In Section 2, we present the SM Lagrangian, and in Section 3 we introduce the effective Lagrangian. 
We discuss Higgs vertices in Section 4 and Z vertices in Section 5. Section 6 gives the relevant partial decay 
widths of the H boson. In Section 7 we list the various H — >• 4f decays. Double Higgs production is discussed 
in Section 8. We discuss perturbative unitarity in Section 9. We give our conclusions in Section 10. 

2 ^ SM : definitions 

In this Section we collect all definitions that are needed to write the SM Lagrangian [24]. The scalar 
field K (with hypercharge 1/2) is defined by 



n>4 (=1 




(2) 



Ref. [23]. 
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The covariant derivative is 

D^K = (d„ - l -g B a ^ a - l -g gl Bl) K (3) 

with g\ = —sq I cq and where % a are Pauli matrices while sq [cq ) is the sine(cosine) of the weak-mixing angle. 
Furthermore 

Wj = ^(*i=F Z^c^-^Bj, A^seBl+ceBl. (4) 

F£v = dn B v-dvBl l +goS ahc B l l l b c v , F® v = d^B®- d v B Q ^. (5) 
Here a, b, ■ ■ ■ = 1 , . . . , 3. The dual tensor is defined by 

f> = ^ Va/5 FV (6) 

Furthermore, for the QCD part we introduce 

Ga -\ a -\ a . rabc b c 

M v = ^gy-dvg^+gsf g^g v . (7) 

Here a, b, ■ ■ ■ = 1 , . . . , 8 and the / are the SU (3) structure constants. Finally, we introduce fermions, 

*=(i) t ^ =l 2^ f 

and their covariant derivatives 



D^L = ( d Li+8 B '[i T i) Yl> » = 0,...,3 

r a = -^x fl , T° = ~ g2 i t ( 8 ) 

£> M y R = (a M + gSj, *,■) Vr, f fl = 0, (9) 



2 V S4 
with g ; - = -se/ceki and 

^ 2 = i-2e„, x 3 = -2e„, x 4 = -2^. (io) 

The Standard Model Lagrangian is the sum of several terms: 

-2sm = -2™ + -^k + «5?gf + «^FP + «£f (11) 

i.e., Yang-Mills, scalar, gauge-fixing, Faddeev-Popov ghosts and fermions. Furthermore, for a proper treat- 
ment of the neutral sector of the SM, we introduce a new coupling constant g, defined by the relation 

gQ = g(i+g 2 r), (12) 

where F is fixed by the request that the Z — A transition is zero at p 2 = 0, see Ref. [25]. The scalar Lagrangian 
is given by 

Sf K = -{D fl KyD ll K-fi 2 K^K-^X(K^K) 2 . (13) 
We will work in the ^-scheme [25], where 

„>=ft-2V, X=\if* (14, 

Furthermore, we introduce v = \/2M/g. and fix fih order-by-order in perturbation theory by requiring 
< 0IHI0 >=0. 
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Table 1 : A selection of relevant d = 6 operators 



^ K = -i(K t K) 3 

< = # 2 (K + K) (D^K^K 

K = i§ 2 {D^Y XaD^KF^ 
4=g(Kt K -v 2 )F; v F; v 
^=gK\,KF« v *J v 
^ e 2 v =g(KtK-v 2 ) Fj y Fj, 
^g=^(K t K-v 2 ) G£ V G£, 
^ f 2 = g 2 (K 1 " K - v 2 ) v L K c b R + h. c. 

< = V L D M b R D M K c + h. c. 



4 = 



^ (KtK) d M (KtK) 



(K^K 



(D M K) + K 



01 = ig 2 (D M K) f D M K^ v 
^ 2 =g(KtK-v 2 )F° v F« v 

^ v =^(Kt K -v 2 ) f; v f« v 

ff\ = g 2 (K^K-v 2 ) \j/ L Kt R + h. c. 
G\ =^r L D M t R D M K + h. c. 



Table 2: Alternative single-fermionic-current a! = 6 operators 



^ f 5 = (K^K) VL^VL + h. c. 
^ = (K t D M K) v R 7^ R + h. c. 
0} = (Kh a D^K) v L x' ! 7^v L + h. c. 



3 Simplified effective Lagrangian 

Our minimal list of d = 6 operators is based on the work of Refs. [26,27] and of Refs. [28-31] (see also 
Refs. [32,33], Ref. [34], Refs. [35-39], Refs. [40-42] and Ref. [43]) and is given in Table 1. 

j2f = jgf SM + £^^f= 6 , (15) 

The structure of the d = 6 operators is chosen in such a way that, with /3/j = 0, no term proportional 
to \/g will appear in the Lagrangian (a part from irrelevant constant terms). Operators containing F^ v are 
CP-odd, the remaining ones are CP-even. 

Additional operators not included in Table 1 have been considered in Eq. (4) of Ref. [38] and are given in 
Table 2. In certain models their effect can be comparable to the one of However, they do not contribute 
to the Hqq vertex, as explained in Ref. [38], because the vector current qy^q is conserved. For a complete 
list of d = 6 operators (other than the four-fermion ones) we refer to Table 2 of Ref. [27]. 

For the single-fermionic-current operators we have adopted the (simplified) choice of Ref. [30], discard- 
ing the chromomagnetic dipole moment operator, which affects the process gg — > ft; in general it is known 
how to remove derivatives acting on the spinors using integration-by-parts. For a complete classification we 
refer again to Table 2 of Ref. [27] where there are 13 operators of dimension six involving single-currents 
of quark fields. 

If one restricts the analysis to the calculation of on-shell matrix elements then there are linear com- 
binations of operators that vanish by the Equations-Of-Motion (EOM). Under this assumption there are 
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redundant operators, e.g. G^, which can be expressed in terms of a d = 4 operator (K^K) 2 , of G K ,Gg K 

1 2 

and of higher dimensional Yukawa interactions involving \j/Y an d three ^-fields, i.e. G f ' . Since we are 
working with unstable particles, the use of EOM should be taken with due caution; indeed, only S -matrix 
elements will be the same for equivalent operators but not the Green's functions. 

It has been pointed out in Ref. [44] that, even if the S -matrix elements cannot distinguish between two 
equivalent operators G and G' , there is a large quantitative difference whether the underlying theory can 
generate G' or not. It is equally reasonable not to eliminate redundant operators and, eventually, exploit 
redundancy to check S -matrix elements. If one eliminates them, whenever the Higgs boson is taken on-shell 
and the full set of d = 6 operators of Ref. [27] is used, the presence of G^ is redundant, and one should set 
4 = 0. Strictly speaking, the last statement only applies to single-Higgs processes; the argument is simple 
(see Appendix. D of Ref. [44]), given a theory with a Lagrangian Jzf [<j>] consider an effective Lagrangian 
,if eff = + g & + g ' G' where G- G' = F[ty] 8^/8^, and F is some local functional of <|>. The effect of 
G' on Jzfgff = S£ ' + gG is to shift g — > g + g' and to replace § — >■ <|) + g 1 F and F contains terms with several 
fields, Q.E.D. 

The effective Lagrangian of Eq.(15) is one possible way of parametrizing deviations of the Higgs cou- 
plings to SM particles; if confirmed, these deviations require new physics models that are the ultraviolet 
completion of the set of d = 6 operators. However, there are specific assumptions in considering Eq.(15), 
namely decoupling of heavy degrees of freedom is assumed and absence of mass-mixing of the new heavy 
scalars with the SM Higgs doublet. 

We postpone a more detailed discussion of non-decoupling effects to Section 6.5; here we note that 
Eq.(15) comprises all heavy physics effects at scales below A, and in a decoupling scenario A is the mass 
of the additional, heavy, degree of freedom. A typical non-decoupling scenario is given by the inclusion of 
a scalar triplet; here higher dimensional operators are not suppressed by inverse powers of the triplet mass. 
It is considerably more difficult to construct a perfectly sensible low-energy effective theory in the non- 
decoupling scenario and the construction is model dependent, e.g. it has been shown in Ref. [45] that (in the 
heavy triplet case) A is related to the ratio of the renormalized triplet VEV to the renormalized doublet VEV. 
Therefore, additional work is needed in handling models showing a non-decoupling behavior, e.g. looking 
for the presence of alternative large parameters. 

3.1 From the Lagrangian to the S -matrix 

There are several technical points that deserve a careful comment when we construct S -matrix elements 
from the Lagrangian of Eq.(15). 

• Field-scaling, parameter re-definition 

We perform field re-definitions so that all kinetic and mass terms in the Lagrangian of Eq.(15) have the 
canonical normalization. First we define 



H = H 

M 2 



1 + ^ (4 + 4 + 2^ K ) 



*° = * [ l + -g; (4+4)], ^ = ^[1 + ^-4 

then we introduce new parameters, 



(16) 

M 2 



(17) 



M = M ( 1 + ^2-4)' c °=ce I 1-^2-4)- (18) 
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Finally we rescale again the fields 



A 2 



redefine the weak-mixing angle as 



M 2 
A 2 



Cq=Cq [ 1 -A'-^-SQCgal ) , 



and introduce Higgs parameters 



M 2 = 



M / 1 3 „ \ 



— , M 

M 2 - 16— -yflK, 



JSh= [i-^T (24 + 4 + 2a aK )]j3 H - 4^ a K . 



It is worth noting that a different definition of G\, i.e. 

^ =g 3 K t K^ L Kt R + h. c., G\ = g 3 K\j/ L K c b R + h. c, 

requires a re-definition of the t — b bare masses, 

— 2 — 2 

M t =M t -2V2-^aJ, M b =M b -2V2^ r ^. 

In the option of Eq.(l) the Hff-Yukawa couplings are 

1 M 2 



&BSL ~ -^Hff ~ 2 A 2 " 



1 M 2 r M\, 1 3 

2 A 2 " ( a K + 2 ^K 



(4 + 2a 3K ) -4>/2aj 
Hbb 



Htt 



(4 + 2^ K )-4^2 £ 

while, following Eq.(23) and Eq.(24), we obtain 

1 M t M / , \ __ 1 M b M / , \ __ 

^kff = ^-2«-^- (4 + 2fl5 K )Htt--g-^- (4 + 2^k) Hbb. 

• gauge-fixing term 

We define a modified gauge-fixing term for the W,Z -fields, 

c ± = -a M wJ+^ w M^ ± , c z = -d M z M +^ z ^V 



where the gauge-parameters are 



It is straightforward to show that 



1 



E M. 



C ± = -^Wj + Mc^±, c z = -^d M Z M +| z ^o, 
where the gauge parameter is 

M 2 _ , 

Note taht the photon gauge-fixing term remains unchanged, i.e. 

c A = -a M A 



(29) 



(30) 



(31) 



With our choice for the scaling factors, the parameter redefinition and the form of the gauge-fixing term it 
follows that the part of the Lagrangian which is quadratic in the (bosonic) fields reads: 

J^ 2 bos = -^W+^W- -M 2 W;W M - ld M Z y d M Z v - \ ^Z^Z^ - I^A V ^A V 





+ 4^ c e s e 



2 H FT FT T T 2 FT FT 2 

4 (a M Z v 5 V Z M - ^A v 5 V A M ) - 2e' lv ^ « e 3 v (^Z v d a Z„ - ^A v d a A^ 

+ 4(1-2^)^ [4 (^Z v ^A v - a M Z v 5 V A M ) + 26^ a e 3 v ^Z v d„A J (32) 

Therefore, kinetic and mass terms are SM-like, and the bare p parameter is ff (l/A 2 ), by construction. 

• Dyson resummed propagators 

Dyson resummed propagators are crucial for discussing several issues, from renormalization to Ward- 
Slavnov-Taylor (WST) identities [46^48]. Consider the W or Z self-energy; in the SM we have 



<v (p 2 ) = n v v (p l ) 8, v + n; v (p')p, Pv . 

Once d = 6 operators are added the W Dyson resummed propagator remains unchanged, i.e. 



t vv ^„2n 



t vv ^„2x 



—WW 

V 



■ + ■ 



„ww 



^ 2 +m 2 -< w (V + m 2 -< w ) (^ 2 +M 2 -< w -nr^ 2 ) 

while the Z propagator changes as follows: 



(33) 



(34) 



For the if propagators we get 



IT 



V 



zz 



T zz , „ M 2 9 e 3 



2 , T 2 M 2 ' 



l 



/? 2 +M 



2' 



(35) 



(36) 



with the gauge parameter defined in Eq.(30). 
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T"_ M T"_ M 

'Z 'ZZ^ 




Figure 1 : Example of Ward-Salvnov-Taylor identity; the grey circle denotes insertion of d = 6 operators, black circles 
denote the replacement of the polarization vector by i times the momentum flowing inwards. Z and <j)° lines represent 
Dyson resummed propagators. 



• WST identities 

With the Feynman rules developed above we can prove WST identities; we show an example in Figure 1 
where one should take into account that all lines must be on-shell otherwise there are additional terms 
involving FP-ghost lines, i.e. BRST-invariance requires also effective operators involving ghost-fields. 

• Wave-function factors 

Due to the rescaling of the fields each external leg in a S -matrix element has to be multiplied by a factor; 
the argument is general, given a Lagrangian 

Jz? = Z~ 2 (|)A~ 1 (t> + /(|) (37) 

we have to normalize the source / in such a way that the residue of the two-point S -matrix element is one; 
therefore we fix J — > Z~ l J for the S -matrix element containing one external (f> -line. We define Z, = 1 + 5Z, 
and obtain 

<5Z H = ( a K + a K + 2a dKj , 5Z 4>° = -TT ( a K + a l) ! 5Z *=^r a L ( 38 ) 

a s 7 .M 2 ^ 3 M 2 Sg 3 

<5Z W = 0, 8Z z = -4— T — ai / , dZ A =A— T — a 3 w . (39) 
A 2 c e A 2 co 

3.2 Nature of d = 6 operators 

The d = 6 operators are supposed to arise from a local Lagrangian, containing heavy degrees of freedom, 
once the latter are integrated out. Of course, the correspondence Lagrangians — > effective operators is not 
bijective since many different Lagrangians can give raise to the same operator. Nevertheless these operators 
are of two different origins [49] : 

• T -operators are those that arise from the tree-level exchange of some heavy degree of freedom 

• L -operators are those that arise from loops of heavy degrees of freedom. 

The L -operators are usually not included in the analysis. The accuracy at which results should be presented 
is given by 

where LO means the first order in perturbation theory where the amplitude receive a contribution. To be 
precise, if the underlying theory is weakly-coupled operators containing field-strength tensors cannot be 
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T -operators, and their Wilson coefficients are 1/(16 n 2 ) suppressed. If only T -operators (coefficients of 
& (1)) are included only 14 out of 34 entries in Table 2 of Ref. [27] are relevant. There is another caveat: 
d = 6, L -operators have Wilson coefficients ~ 1/(16 7T 2 ) and d = 8, T -operators are ~ v 2 /A 2 ; therefore, 
below 3 TeV one should include both of them or none of them. 

3.2.1 Insertion of d = 6 operators in loops 

The question remains on insertion of d = 6 operators in SM loop diagrams. This is better discussed in 
terms of a concrete example: consider the Lagrangian [29] 

Se = ifsM - \ d^Sd^S - -Ml S 2 + Ms K 1 " KS, (41) 
where S is a heavy (scalar) singlet The interaction is 

JSJU = ^ Ms (H 2 + <S?%° + 2 V) 5. (42) 

In the limit M s ->»we have 

^^s L M + ^KtK+^ K . (43) 

The <i = 4 operator in Eq.(43) can be absorbed through a parameter redefinition, and we are left with a 
contribution to the d = 6 operator dK . Clearly, & dK (as well as &^ and ^|) is a T -operator [29,50,49]. 



H 



P2 



'■ " *•• <- Pi 

Figure 2: The three-point function H 3 with the insertion of the 0j K operator (left) and the same conttibution in the 
full Lagrangian of Eq.(41). 



Imagine we want to compute the H 3 Green function: we analyze the ultraviolet (UV) behavior of the 
two diagrams in Figure 2. In the effective theory (left diagram) there is an UV divergence and one option 
would be to subtract it by introducing counterterms in ££d=(>- However, this shows how the insertion of 
local operators of higher dimensionality in SM diagrams is not really consistent since, in the full theory, the 
corresponding diagram is not divergent. If we introduce A 2 = the diagram behaves like A~ 2 In A, 

i.e. the divergence is controlled by the heavy mass. From this point of view it is important to stress that 
one should avoid using a cutoff procedure with the dimensionful parameter A. Computing with dimensional 
regularization gives a different pole structure reflecting the different counterterms in the full and effective 
theory. This difference is independent of infrared (IR) physics, since both theories have the same IR behav- 
ior. As we have seen, there can also be logarithmic dependence on A; if these logarithms are included they 
must be summed. 
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To give an example we consider again the two diagrams of Figure 2 with s = — (p\ + pi) 2 . In the 
effective theory the insertion of 0g K (left diagram in Figure 2) produces 



4ff : 



3 g 



{q+p\ 



(44) 



where n = 4 — e. Suppose that we use a cut-off regularization, the integral is G{\) for A — > oo but the same 
is true for all integrals containing the insertion of <^ K operators; therefore, all these diagrams are of the 
same order and cannot be neglected. In dimensional regularization (DR) we obtain 



,dr 

'eff 



3 

4 g MA 2 



-s-3M 2 



1 , Mr 
- -ln^- 

£ s 



+ finite part 



(45) 



where l/e = 2/(4 — n) — y — \nn and /Xr is the renormalization scale. In principle, we could add countert- 
erms (in the MS scheme) to remove the UV pole and make a choice for the scale, /^r, which minimizes the 
remaining logarithms in the UV finite part. After subtracting the UV pole we can say that the insertion of a 
d = 6 operator produces a result 



rren 



d=6 



M H 



The insertion of a d = 8 operator, always working in dimensional regularization, gives 



M H 



(46) 



(47) 



etc. Note that, with cutoff regularization, both integrals would be of <^(1). Note that in a mass-independent 
scheme like MS the conditions for the decoupling theorem are not satisfied. Furthermore, the logarithms of 
the renornalization mass may become large. In principle, the problem can be solved but the solution requires 
matching conditions (for a discussion see Ref. [51]). 
With the full theory at our disposal we compute 



3 M2ji| 



/fun = 



M 



i 



Working (for simplicity) with M \ <C s <C Af| we obtain 



(48) 



3 M^jif 



/fun = ^ 



Ms 



C(2)-Li 2 1 



s + iO 



(49) 



We can identify A = M 2 /n$, expand in s/M 2 , and obtain 



/full = ~ V 



[1--- 

MM 2 4M 2 




(50) 



The first term in 7f u n reproduces the d = 4 operator of Eq.(43) while the second term corresponds to the 
d = 6, {?g K operator. There is no UV divergence in 7f u n and the logarithm is uniquely fixed. 
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An alternative way to understand the two different approaches is the following: we start from Eq.(48) 
and expand in the integrand 



1 



1 



1 (g+Pi) 



Ml 



{q + pi f + M 2 Ml 
which is equivalent to inserting d > 4 operators or introduce Feynman parameters: 

J = J d n c 



(51) 



1 



q 2 ((q+Pif+Ml) (q+pi+pi) 

[ dx [ dy Ml (x — y) -sy (1 —x) 
Jo Jo L 



(52) 



use a Mellin-Barnes representation, and expand as follows (Ms — > °°): 



dv [Ml 



1 

2ni J- 
2m Ml 



v-1 



dx I dyB (v , 1 — v) (x — y) v 1 y v (1 — x 



r 2 ( s )r 2 (i- s ) 



-Mf 



(53) 



Here B(x,y) is the Euler beta-function. Using the well know Laurent and Taylor expansions of the Euler 
gamma-function we obtain the result summing over the poles at s = —n: 



1 



\n 



Ml 



n+l 



1 



Ln + 1 



(54) 



The result is manifestly UV finite, term-by-term, and has the correct structure of logarithms. 

3.2.2 Admissible operators 

Missing a candidate for the BSM Lagrangian, we will not deal with renormalization of composite oper- 
ators; therefore, we will not include local operators in loops. To be more precise we will use the following 
set of rules: 

1. operators altering the UV power-counting of a SM diagram are non-admissible 

2. operators that do not change the UV power-counting are admissible only in a very specific case: we 
say that a set of SM diagrams is UV-scalable w.r.t. a combination of d = 6 operators if 

• their sum is UV finite 

• all diagrams in the set are scaled by the same combination of d = 6 operators. 

To explain with one specific example, let us consider the HWW vertex with off-shell lines and no wave- 
function factor inserted: 



'HWW 



-gM 



l + (4-24 + 2^ K ) 



M 2 i a M 

A z 



+ 8 4 M r MV _ 4 M p v + 2p v ip £ + p , pl) + l6a l Y ™ £ *^ VpiaP2p (55) 



A 2 



A 2 
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with T^ v = P2 p\ — pi ■ P2 S^ v and P = p\ + p2- Consider the one-loop diagram contributing to H — > yy 
containing a W loop: the operators ^ and (9^ , ^ e ' v change the UV power-counting of the original SM 
diagram and are non-admissible. 

In the one-loop (bosonic) amplitude for H — > yy there are three different contribution, a W-loop, a 
charged if -loop and a mixed W — if loop. We find 



^HW* = ig 1 + («K - 2 «K + 2 «5k) 



+ ^ 



l + (a^ + 2a dK J 



1 + 



4 



M 2 
A 2 " 

,1 



M 2 
A 2 " 



P2 V 



v , . 4 P2-P2 v 
Pi + la K—^~Pl 

■ 4 Pi -P2 v 



A 2 



M' 



2a K + 2a dK ) +^Ma K 



1 Pi -Pi+Pi-Pi 



A 2 



+ gM (a\-a^ + 2a dK 



PP 



-M 2 



A 2 



It is straightforward to conclude that the SM one-loop, bosonic, amplitude for H 
line is UV-scalable w.r.t. the combination 



Cbos — -^2 {a\-2a\ i + 2a dK s j , 



(56) 

YY with on-shell Higgs 
(57) 



which could be admissible. However, in the one-loop amplitude we also have FP-ghost loops with vertices 
(see Eq.(201)) 



1 



M- 



v hx±x ± = -2# M l l + [ a K + 2a dK ) -jp 
Therefore the bosonic component is only UV-scalable w.r.t. the combination 

M 2 



. M / 3 \ 

C bos - -JJ y a K + 2a dYL ) ■ 



(58) 



(59) 



Similarly, we consider the yWW, yWif, y# and yX vertices, which also appear in the one-loop bosonic 
amplitude for H — > yy, and conclude that the latter is UV-scalable w.r.t. the combination 



r 2 

'-bos 



M 2 ^ 
A 2 ? e 



(60) 



which is also admissible. Obviously, the wave-function factors of Eqs.(38)-(39) are also admissible. To be 
more precise, the one-loop bosonic amplitude for H — > yy is made of three different families of diagrams, 
shown in Figure 3. We find that the yyWW, yyWif and YYM vertices are all UV-scalable w.r.t. 2C 2 os . 
Furthermore, the vertex yHW<|) is UV-scalable w.r.t. C^ os +C 2 0S . The underlying algebra is such that the 
quadrilinear vertex with two ys is equivalent to the square of the trilinear vertex with one y (to & (l/A 2 )) 
and the quadrilinear vertex with one H is equivalent (to the same order) to the product of the two trilinear 
vertices, with a y and with a H. As a consequence, there is a non-trivial scaling factor which is admissible, 
not spoiling the UV behavior. 

The fermionic amplitude for H — > yy contains a top-quark loop and a bottom-quark loop. The top 
contribution is UV-scalable w.r.t. the combination 



w. 



fer 



M 2 

1 + [ai + 2a dK ) -^2" 



+ 



1 M 2 



4V2 A 2 



2a\ + 



PP 



M 



2M 2 3 
5 a i 



(61) 
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(W\A/ 




H 



W/(|)/X 



± 




W/(|) 




Figure 3: The three families of diagrams contributing to the bosonic amplitude for H — > yy; W/(|> denotes a W-line or 
a 4>-line. denotes a FP-ghost line 



while for the bottom-quark we have 

b 1 M b r / - \ M 2 ] 1 M 2 r 2 P-P-lMl -, 

Cfer = "2 ^ M I 1 + K + 2 ^kJ AT] + ^ aT [ 2 «f + m 2 af \ ■ (62) 

One example of L -operator is given by in Figure 4 with contributions from heavy colored scalar fields 
transforming in a (C, T , F) representation of S£/ (3) <g> SI/ (2) (g) £/(l), e.g. the (8,2, 1 /2) representation 
of [52-55]. Since the additional colored scalar (weak-isospin) doublet contains also an electrically charged 




Figure 4: Example of diagram giving a contribution to the d = 6 operator of type L. Solid lines represent colored 
scalar fields, e.g. transforming in the (8,2, j) representation of SU(3) <g> SU(2) (g) t/(l). 

scalar (and two neutral scalars) it will contribute to the decay H — >• yy [56]. As long as the scalars are in a 
representation (C ,T ,Y) such that C C B 8 there will also be a contribution to gluon fusion. 

3.3 Effective theory and renormalization 

There are two conceptual frameworks to discuss renormalization and effective theories. In one case we 
are only interested in setting up an expansion in power of E/A where A is the cutoff and E is the scale 
relevant for a given set of processes. 

Counterterms are introduced to remove UV-divergences and, in presence of d > 4 operators, an infinite 
number of them is required. However, once the requested precision of the calculation is fixed only a limited 
number of term is needed. 

This is not the goal for Higgs physics where we want to search for new physics without committing to 
a particular extension of the SM. The effective theory should simply capture the low-energy effects of the 
underlying, BSM, theory and must be replaced by a new one when E is approaching A, where it should be 
discarded. Having this difference in mind we proceed in discussing renormalization. 

The processes gg — > H and H — > yy are special in the sense that there is no tree-level coupling and 
therefore the NLO (one-loop) amplitude is UV finite. This is not the case for other processes, i.e. H — > bb 
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etc. In general we will have 



^ = /(Kc}) [A L o({po})+A mjo ({po})\+A AC ({a AC ,p }), (63) 

where {po} is the set of bare parameters (masses and couplings), {a AC } a set of effective parameters; fur- 
thermore A L0 (A NL0 ) is the LO(NLO) SM amplitude. Since Anlo contains UV divergences we introduce 
counterterms 

PO = Pren + SZ p , (64) 

where p ren is the renormalized parameter and 8Z p contains counterterms. If A' denotes the derivative of the 
amplitude w.r.t. parameters we obtain 

^ = f({a A c}) A LO ({p ISD })+A' LO ({p ien }) ® {Z p }+A NLO ({p ien }) +A AC ({a AC , j p re n})- (65) 
The combination 

Alo ({Pren}) ® {Z p } + A NL0 ({pren}) (66) 

is now UV finite. Note that we have replaced po — > p Kn in A AC because in the full theory A AC is of the same 
order of A NL0 , i.e. renormalization of a AC can only be discussed in the context of the full theory. In a sense, 
the a AC parameters are already the renormalized ones. 

There is a final step in the procedure, finite renormalization, where we have to relate renormalized 
parameters to physical quantities (e.g. e 2 = g 2 s^ = a /(An)), 

Pren = Pexp + F ({Pexp}) • (67) 

This substitution induces another shift in the amplitude 

A LO {{Pren) -> A LO ({Pexp) + A L0 ({/?exp) F ({/?exp}) , (68) 

with p ren = p exp in both A NL0 and A AC . This set of replacements completely defines our renormalization 
procedure. 

A subtle point is the following: in the process H — > yy we have a bosonic component of A L0 and a 
fermionic one and both are UV finite. Therefore, as long as all tree-vertices in the bosonic part are scaled 
with the same factor, we would like to have 

* = /bos ({a AC }) A L b ° s ({/? }) +/fer (Kc}) A% ({p Q }) +A AC ({a AC , Po }) . (69) 

LO implies one-loop diagrams where the splitting bosonic-fermionic has a meaning. Once we try to go to 
NLO (i.e. two-loops) the splitting is not definable and renormalization is requested, i.e. one has to insert 
counterterms in the one-loop diagrams. Clearly, an arbitrary scaling of the two LO components kills two- 
loop UV finiteness (at least in the electroweak sector). The two-loop electroweak corrections to H — > yy 
are —1.65% at Mh = 125 GeV [57], therefore neglecting them is tolerable but the internal inconsistency 
remains. The effect on gg — > H is larger, 6 (5%). 

To conclude this section we compare the BSM scenario with heavy degrees of freedom and the SM one 
in the limit of infinitely massless top-quark. In this case we have a coupling Hgg of the form 

-£/ = --A, SM HG£ V G^ V , (70) 

where A, SM has inverse mass dimension. The important point is that A. SM is computed by matching the effec- 
tive theory to the full SM [58-69]. Even more important is the fact that A. SM in the effective theory is the 
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renormalized one, with its renormalization constant computed to all orders [70]. Therefore, the logical steps 
are: first renormalization in the full theory, then construction of the effective one. 

To be more precise we consider a theory with both light and heavy particles; the Lagrangian is j£f (ra) 
where m is the mass of the heavy degree of freedom. Next, we introduce the corresponding Jz? e ff , the effective 
theory valid up to a scale A = ra. We renormalize the two theories, say in the MS -scheme (taking care that 
loop-integration and heavy limit are operations that do not commute), and impose matching conditions 
among renormalized "light" 1PI Green's functions 



For the case where the full theory is the SM and ra = Mh the whole procedure has been developed in 
Ref. [71]. 

4 Higgs vertices 

We are now in the position of writing the complete expression for vertices. There are different level 
of implementation and accuracy. We start with LO-inspired accuracy where the SM vertices are at LO and 
the tensor structure of the vertices is the same as the LO SM one but every coefficient coming from the 
effective Lagrangian is kept. Next we go to LO-improved accuracy where extra tensor structures from the 
effective Lagrangian is included. Finally there is an NLO-inspired accuracy where the SM components are 
at NLO but contributions from d = 6 operators are included only under the constraint that they do not spoil 
UV-finiteness. With the introduction of the following tensors we obtain 



iftdi(M) = r») 



jU < m. 



(71) 



afijiv 



(72) 



P 






(73) 



P 




HZZ 







(74) 
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Vp 2 



HAZ 



[1-2Sq) a w + 8c 2 e $e (ay -ay 



( " 5 , » 4 

I cga K + sea K 



M 

7^r e 



n M 

M 



g (l - 2f|) a^y + 2c e *e (a^v - aj w 



A 2 c e 



a 5 i a 4 \ M v 



(75) 




Vp 2 



H W W 



-gM 



1 



M 2 

2a^-a^-2ag K ) — j 



8^ 



(76) 



p 


/ 










^vbp2 


Hgg 


M 

A 2 e 


(77) 




Vp 2 



Hit 



Af t 
M 



1 + (aK + 2a 5 



M 2 ' 
k) A 2 " 



„ /-M 2 j 1 P 2 3 

A 2 4srt A 2 1 



(78) 



Similarly for f = b we have 



Mb r 

M 

2 



„ /rM 2 1 P 2 4 

A 2 4\/2A 2 



— M 2 

H(P)^b(pO+b(p 2 ) = -^^ l + (a| + 2a 5K ) 



(79) 



16 



5 Z couplings 

The Zff vertex can be parametrized as follows: 

hf(l + f)-2Q f K f s 2 e 



2ce 



where % is the third component of isospin and Qi = — 1, Q v = 0, Q u = 2/3 and <2d 
part reads as follows; 



(80) 



-1/3. The anomalous 



Ap f 
AfCf 

6 Partial decay widths 



M 2 
_M 2 



A 2 & 



a^-32Q f I 3i (\-So))cl ay 

4+4 (l+4j2f/ 3 f*fl(l-*e)) *Ve4 



(81) 



In this Section we compute the partial decay widths of the Higgs boson for the most relevant channels: 
first we introduces the dimensionless coupling 



1 



■ 0.085736 



TeV 



GpA 2 I A 

which parametrizes deviations from the SM results. Furthermore, we introduce new couplings 

ga Y =A Y , ga\=A\, g 2 a\=A\, ga g =A g 
g 2 a l K =A^, g 2 a\=A\, g 2 a dK =A dK , 



gOf 



1 ■ —A 1 - 



4y/2 M 



2 1 M b a2 



4V2 M 



(82) 

(83) 
(84) 
(85) 



and express all amplitudes in terms of a SM-component (eventually scaled by the effect of d = 6 operators) 
and by a contact component, as shown in Figure 5. We introduce an auxiliary coefficient, 



A 3 



A^ + 2- r +4A 3K - 



(86) 




Figure 5: Amplitude for a two-body decay of the Higgs boson (dash line) including LO+NLO SM contributions with a 
sum over all one-loop diagrams (i); SM diagrams are eventually multiplied by a universal scaling from d = 6 operators 
(black circle); the grey circle represents a contact term. 

We will now show results for various decay processes. 
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• H — > YY 



For H — > YY the SM amplitude reads 



where 



v M 

.Pi Pi 



= F SM \8^ + 2^)e,( Pl )e v (P2) 



Sr w 1 l M 2 

-g M F SM - -g T =-F su --g==-F SM . 
2 M 2 M 



(87) 



(88) 



•<SM 



Ml 



6 + =^- + 6 (M 2 -2M 2 ) C (-M 2 ,0,0;M,M,m) 



FL = -8 - 4 (m 2 - 4M 2 ) C (-M 2 , , ; M t ,M t ,M t ) 
etc. The Co function is given by 



Co (— M 2 , 0,0; m,m,rnj = — 



1 



2M 2 



In 2 



) l ' 2 + l 



(i-x)!/2_r 



(89) 



(90) 



where x = 4 (m 2 — iO) /M 2 . Note that there is no need to split the result for this Co -function into the two 
regions x > 1 and x < 1 since the iO prescription uniquely defines the analytic continuation. We find 



^r H -m = (4V2G F ) 1/2 {-f [c^F^ + S^Cr^ + S^C^^] +F AC } 



g6 

V2 



M 2 (f 2 A^ + 4A 2 +c e feA^) . 



where the scaling factors are given by 



-07 

-w 



for the W -loop and 



C 



YY 



-07 



>1 



1 + 
1 + 



4^2 



8A^c e (fe + - )+A^-A f 



1 

50 



} 



(91) 

(92) 

(93) 
(94) 



for the quark loops. 

The amplitude is the sum of the W, t and b SM components, each scaled by some combination of 
Wilson coefficients, and of a contact term. The latter is & (g$) while the rest of the corrections is 6 g$)) . 
However, one should remember that O l y are operators of L-type, i.e. they arise from loop correction in 
the complete theory. Therefore, the corresponding coefficients are expected to be very small although this 
is only an argument about naturalness without a specific quantitative counterpart (a part from a 1/(16 K 2 ) 
factor from loop integration). 

The result for H — > gg follows straightforwardly. 
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• H^gg 



The result for H — > gg is straightforward. Including also the b -loop we obtain 

1/2 r a s (M 2 ) 



= (4V2G F y 



71 



C gg /4 + CfF s b M )+^M 2 A ? 



where the scaling of the quark components is given by 



Cf g = ^M 2 
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-Mi 



1 + 



1 + 



86 

4^2 

86 

4V2 



(Ai-Al) 
(Ai-Af) 



(95) 

(96) 
(97) 



• H^bb 



For the H — > bb amplitude we have to examine again if the are UV-scalable diagrams. In this case there is 
a tree-level amplitude, and renormalization is required. At NLO there are two type of diagrams, the abelian 
ones involving the Hff vertex and the non-abelian ones involving a HVV (HV^Hcjxj)) vertex. Therefore we 
have to search for the unique combination that multiply all the vertices, which is 



M 2 



The SM amplitude reads as follows: 



ySM _ 
"H^bb 



,3 M * 77 SM 
? ■=" ' , 



M 



H^bb 



u(p 2 )v(pi). 



(98) 



(99) 



The expression for can be found in Section 5.9.4 of Ref. [24]. Renormalization and QCD corrections 

are discussed in Section (11.2—11.4) of Ref. [24]. The complete amplitude for H — > bb is 



^bb = (4^Gp) M b u(p 2 )v( Pl ){^C bb F^- hb 
-M 2 

=|A4-16(A| + 2A aK +A f 2 )]}, 



C 1 



bb 



128 y/2 
1 



1 + 



86 



Ak+A 3 K + 6A dK j 



(100) 



(101) 



In the SM, NLO corrections to the amplitude include a QED part so that, technically speaking, the process 
is H — > bb(y), i.e. real corrections are added. There is also a contribution from the d = 6 operators 
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(102) 



3V2 A 2 ' 

which, however, is not infrared divergent and will not be included. 

In this Section we have considered partial decay widths of the SM Higgs boson; in the SM, the common 
belief is that (for a light Higgs boson) the product of on-shell production cross-section (say in gluon-gluon 
fusion) and branching ratios (zero-width approximation or ZWA) reproduces the correct result to great 
accuracy. The work of Ref. [72] shows the inadequacy of ZWA for a light Higgs boson signal at the level 
of 5%. Therefore, one should always implement the results of this Section within a consistent off-shell 
formulation of the problem. 
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6.1 gg^H 

In ZWA the inclusive cross section for the production of the SM Higgs boson in hadronic collisions can 
be written as 

a(s,M^j = £ [ dx x [ dx 2 fi/ hl (x l ,^)f j / h2 (x 2 ,^) x 

x f'dzS (z- — J zO {Q) Gij (z;aM),Ml/^;Ml/^) , (103) 

JO \ SX\Xl I \ / 

where yfs is the center-of-mass energy and jU F and jIIr stand for factorization and renormalization scales. 

In Eq.(103) the partonic cross section for the sub-process ij — > H + X, with i(j) = g,q,q, has been 
convoluted with the parton densities f a /i % for the colliding hadrons h\ and hi- The Born factor reads 



a (o) = Gf 



/SM 

q 

q=t,b 



where Gp is the Fermi-coupling constant; the amplitude is generalized to 



(104) 



jg= £ c f ^ q SM + ^ AC , (105) 

q=t,b 



where the last term is induced by the operator G % , and where the scaling factors are 

J6( A \ l- A}) cf = \ + ^ 



Since ^# t SM and are separately UV finite it is possible to include NLO(NNLO) QCD corrections even 
in presence of anomalous scaling factors The coefficient functions Gjj of Eq.(103)can be computed in QCD 
through a perturbative expansion in the strong-coupling constant as, 

G tJ (z;aM),M 2 u/ntM 2 u /nl) = a^i) £ <%> (z;M*/ , (107) 

with a scale-independent LO contribution given by 

Gfj ) (z) = 8 ig 8 jg 8(\-z). (108) 

The NLO QCD coefficients have been computed in Ref. [73], keeping the exact M t and Mb dependence. 
NNLO results have been derived in Ref. [74] in the large M t limit (see Ref. [75] for the NLO case); analytical 
expressions can be found in Ref. [64]. The accuracy of these fixed-order computations has been improved 
with soft-gluon resummed calculations [76-78]. 

QCD corrections cannot be implemented in the additive part of Eq.(105). To do that one needs a model 
for G % , as done in Section 2 of Ref. [54] where the SM is extended to included colored scalars, so that one 
has 

^ = £.# q + £..#s, (109) 

q S 

where fermions and scalars transform according to some (C,T ,Y) representation of SU(3) (8) SU(2) ® 
U(l), as long as C <8> C 3 8. Complete QCD corrections for fermion and scalar amplitudes have been 
computed in Ref. [54]. 
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6.2 Simplified scenario 

1 3 

If we restrict the scenario to bosonic T -operators only (ax, and aj K ) the scaling factors are: 

1 



t b 4 



1 + ^A° 



4^2 



1 



1 + ^A? 



C bb = 



4\/2 



1 + 



r gg - — m? 
Cb " 16 Mb 



1 + ^A° 



4^2 



§6 



2y/2 L 4V2 



The contact terms are all zero but 



H-s-bb ~~ 



(a^+A| + 6A 5k ) 
M b (a^ + 2A^k). 



(110) 

(111) 
(112) 

(113) 



In this case it is not possible to differentiate bosonic loops from quark loops. 
6.3 BSM Lagrangians 

By BSM Lagrangians we mean those Lagrangians containing new, heavy degrees of freedom that can 
produce d = 6 operators when the heavy particles are integrated out. One of the most important questions is 
about the sign of the Wilson coefficients a, in Eq.(15), i.e. to find the set of coefficients such that 



(114) 



Before entering the discussion on BSM Lagrangian we recall few, well-know, facts about tree-level custodial 
symmetry. The SM Higgs potential is invariant under SO(4); furthermore, SO(4) ~ SU (2)l <8> SU (2)r and 
the Higgs VEV breaks it down to the diagonal subgroup SU (2)y. It is an approximate symmetry since the 
C/(1)y is a subgroup of SU (2)r and only that subgroup is gauged. Furthermore, the Yukawa interactions are 
only invariant under SU (2)l <8> C/(1)y and not under SU (2)l (8) SU (2)r and therefore not under the custodial 
subgroup. Therefore, if we require a new CP-even scalar, which is also in a custodial representation of the 
group, the W/Z -bosons can only couple to a singlet or a 5-plet, as discussed in Ref. [79]. If (A^l,^Vr) 
denotes a representation of SU (2)l <S> SU (2)r, the usual Higgs doublet scalar is a (2 , 2), while the (3,3) = 
1 © 3 © 5 contains the Higgs-Kibble ghosts (the 3), a real triplet (with Y = 2) and a complex triplet (with 
Y = 0). The Georgi - Machaceck model [80] has EWSB from both a (2 , 2) and a (3 , 3). 
To introduce the discussion on BSM Lagrangians we define the following quantity: 



AC = g 6 A 



K- 



(115) 



Assuming A v = and requiring that the coupling HW + W in the decay H — > yy has the standard value, i.e. 
that 

(116) 



*-w — ^w 



we obtain the condition AC = 0. We now examine different models and explicitly compute the corresponding 
value for AC. At the same time we address the question of models allowing for non-standard coupling Htt 
in the loop for H — > yy. 

In general, the basis for a representation of SU (2) can be characterized [81] as a tensor field 



11 — („ 



Jnl 



(117) 
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where G are 51/ (2) -matrices. An irreducible representation of spin n/2 is characterized by a totally sym- 
metric field with n indices. The hermitian conjugate \|/J _ ; . transforms according to the complex conjugate 
representation, \\r' 1 and indices can be lowered using the metric tensor To define the covariant deriva- 
tive we introduce 



J«l l n 

i'l — ('„ 



r=l 






°;;:;;:;;=r1^(-r«)" n 



r=l 



IT 



1] ... (; ...!„ 



n 



The covariant derivative is 

(7) 1 ,y) /, -: / ' = {if 1 V) ' ' ' h, d u +g Wf. I 7 ', 1 " 7 7 f //' •••'}•••'" 
v A I> f/fi...f„ \ n-h )[■■■]'„ M n j[-f n JLi H ...i,...i n 



1=1 



fl 

" w /-'r /, 2- ( - /.../:.../. .../;. . 



where a = 0, . . . , 3, W^' 2 ' 3 = fij; 2 ' 3 and W° = gl fl° . 
Here are a few examples of BSM Lagrangians. 



(118) 
(119) 



(120) 



• Example 1 

Consider the following Lagrangian [29] : 



(121) 



^ = -^^S-iMlS^MsK+KS- \ {D^) a {D^) a - l -M 2 ^\ a 



+ MxK t x a KT 1 a - (D^) ta (D^) a -iiff^V+ ^K^KT + h-c 



(122) 



where S is a scalar singlet and r|,£, are scalar triplets with different hypercharge, see Refs. [16,82]. To be 
more precise, n can be written as a complex symmetric tensor of rank two and £ as a traceless tensor. In our 
case we introduce 



Mn 2 



G F M4' 



GfM^ 



4' 



Projecting onto the d = 6 operators we obtain 



AC 



The scenario of Eq.(116) has a solution 



X s = 



which requires the condition 



s2 



1 I X T] -2X s -2 | 1 



«2 



1 1 



?2 



?2 



X 



■5' 



X,>2 



2 + f 2 



2^" 



(123) 



(124) 



(125) 



(126) 
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• Example 2 

Alternatively, we could consider a Lagrangian [29] 

=£-2 = -S-SM + -S-v 



(127) 



1 



1 



= ~ 4 V V MV - -M 2 V M V^ - ig v V M [((D M K) T K - K^D M K 



which contains 7 = and 7=1 new vector fields; introducing 



gy,u 
G F M 2 V 



we obtain that the scenario of Eq.(116) requires 



~2 



x u _ 8-^-PX v . 



(128) 
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• Example 3 

A mixture of vector and scalar fields [29], e.g. 

=£?3 = =^SM + -^sx 



(131) 



gives 



= - \ dpSduS - ^M 2 S 2 + Ms K^KS - Mvs d M S, 



1 



- -v„„v 

4 
1 



- - U£ v U% v - - M 2 V£ V£ -Lgu V" [ K) + x a K - K + x fl D M K 



AC = ix TI -4|x v + 4X s fl " 



2^ u " « 



M 2 



(132) 
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The scenario of Eq.(l 16) requires large values for /Xys • When we include all scalar an vector fields, Eq.(l 16) 
is satisfied by 



X u= 4 



ci 



2^X V + X n -6X 6 -2X S 1- 



L 5" 



Mvs 

a4 



(134) 



Example 4 



1 2 

In order to differentiate the bosonic amplitude from the fermionic one we need & i ' . One way to introduce 
them is to consider an additional Lagrangian, 



«S?4 — =^sm + 



(135) 
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where % is a doublet 



&z = ~\ {D^D ll %- l -M 2 1 xh+ k (K ! K) (K+x) +h.c. 



+ 



Y x Vl X c t R + y x V L X b R + h.c. 



(136) 



which would produce of the order of (Y X X X )/M 2 . 

Finally, we examine the possibility of a non-zero F AC in Eq.(91). This requires Gy operators. One option 
is to include colored scalar fields [56] but we could also include a real triplet [80,82,83] 



with hypercharge Y = 0. The Lagrangian reads as follows 

JS? 6 = - (D^D^-Mj^ + ^ (tfK) (^) , 
with co variant derivative = — igB^ T a and 



(137) 



(138) 



1 



'0 1 0' 



\ / 2\ 10 / 
which gives the following couplings: 



r 2 = -L|, <>-,■ 



'0 -i 
i -i 
i 



T 3 



10 
00 
00-1 



(139) 



and produces a loop of "% scalars in the H77 coupling. 

Additional examples of BSM Lagrangians can be found in Refs. [84,85] and in Ref. [86]. General 
studies can also be found in Refs. [87,88,8]. 

6.4 MSSM 

In this paper we assume that the starting point in comparing theory with data is the SM. Another choice 
could be to start from the Minimal-Supersymmetric Standard Model (MSSM); in this case all the amplitudes 
should be replaced, e.g. 



(h -»• 77) 

= JZ SM (h -»• 77) +g hH + H - 775- A o ( t h j 



+ £iv/ e ^ hf Aq ^ j +£ g + A, (x Zi ) 



(140) 



withx ; =M2/(4Mf) and 



Ai(x) = 3 t + (t-1)/(t) , A (x) = --5 x-/(x 



and 



/(T) = "4 ln Tr^T 



(141) 



(142) 
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Here ghxx is the coupling of h to X = {H 1 * 1 , f , % { }. 

Given the number of free parameters in the MSSM that are relevant for Higgs phenomenology, the 
present experimental information will clearly not be sufficient to fit the MSSM parameters and a further set 
of Wilson coefficients for d = 6 operators. 

An alternative option would be to integrate out the heavy MSSM Higgses (since Buchmuller - Wyler 
basis only has a single Higgs field). By squaring the corresponding MSSM interaction Lagrangian and 
contracting the propagators in all possible ways the coefficients will be calculable. 

6.5 Decoupling 

In this Section we study the problem of decoupling of high degrees of freedom by considering again the 
decay H — > yy. To be fully general we assume the existence of heavy fermions and scalar that transform 
according to generic Rf and R§ representations of SU (3) [56]. The BSM amplitude is based on couplings 

1 Af f , Ms 

Hff =2^SV HSSr =^ <143) 

where the Xs are numerical coefficients (model dependent) and ;Us has the dimension of a mass. The HS + S 
vertex follows from the following choice of the potential: 

V = V SM + 2 (Ml - X s TrS j S +g 2 X s -± (K+K) TrS+S + • • • (144) 

where S = S a T a (T a are the generators in the R$ representation) and where the trace is over color and SU (2) 
indices of the field S, 

e.g. S in the (8 , 2, ±) of SU(3) (g> SU(2) ® U(l). The amplitude reads as follows: 



( H -> YY ) = N{ Xf Q\ Jit + N c s X s (145) 

M 1 



where Q is the electric charge of the particle and N c is the color factor. In the SM we have 

Xf = \, X s =0, N{ = 3 and Rf = 3. (146) 

The amplitudes are 



2 



" x 2 

with %i = M 2 i /(4Mf). In the limit M t -f«we have 



1 



x f + (x f - 1) / (x f ) Jt% = — j Xs ~ f ( Ts ) 



(147) 



/(xO = x,- + | + ^(x 2 ). (148) 



The limit x, — > gives 

.2 



^bsm (H -> YY) -> \ N! X f Q 2 + i Nl X s Q 2 S ^ (149) 



25 



showing decoupling for PS. As stated in Ref. [56] there is decoupling in the theory when v = \flMjg <C 
Ms ; therefore, colored scalars disappear from the low energy physics as their mass increases (Appelquist- 
Carazzone "decoupling theorem" [89]). However, the same is not true for fermions, as shown in Eq.(149). 
We repeat here the argument of Ref. [90] : for a given amplitude involving a massive degree of freedom 
(with mass m), in the limit m->«we will distinguish decoupling A ~ 1/m 2 (or more), screening A — > 
constant (or lnm 2 ) and enhancement A^m 1 (or more). Any Feynman diagram contributing to the process 
has dimension one; however, the total amplitude must be proportional to T^ v = p% p\ — p\ • pi because 
of gauge invariance. For any fermion f the Yukawa coupling is proportional to mf/Mw and T has dimension 
two; therefore, the asymptotic behavior of any diagram must be proportional to r/m t - when nif — > oo. The 
part of the diagram, which is not proportional to T, will cancel in the total because of gauge invariance (all 
higher powers of nif will go away and this explains the presence of huge cancellations in the total amplitude). 
At LO there is only one Yukawa coupling as in NLO(NNLO) QCD where one add only gluon lines, so there 
is screening. 

It is worth noting, once again, that electroweak NLO corrections change the scenario: there are diagrams 
with three Yukawa couplings, therefore giving the net m 2 behavior predicted in [91], so there is enhancement 
and, at two-loop level, it goes at most with m\. At the moment, the NLO electroweak corrections for heavy 
scalar are missing and no conclusion can be drawn on decoupling at NLO. 

In conclusion the decoupling theorem [89] holds in theories where masses and couplings are indepen- 
dent. In all theories where masses are generated by spontaneous symmetry breaking the theorem does not 
hold in general. Another typical example is given by the inclusion of a Higgs triplet: if the triplet develops 
a vacuum expectation value (v^ ) then the p -parameter deviates from unity at the tree-level [92,93] with 

p L0 = 1-2\/2Gfv|, for Y = l 

p L0 = 1+2^2 G F v 2 , for Y = 0. (150) 

We will not discuss details of renormalization but one should always remember that whenever p / 1 at 
tree-level quadratic power-like contribution to Ap are absorbed by renormalization of the new parameters of 
the model and p is not a measure of the custodial symmetry breaking [94]. Alternatively we could impose 
custodial symmetry, = v n , in a model with both triplets; an example is found in Ref. [80] containing 
SU(2) L ® SU(2) R multiplets. 

As far as the triplet contribution to H — > yy is concerned it is also known [94] that decoupling occurs 
only for special values of the mixing angles in the triplet sector. 

An important tool in studying decoupling of heavy degrees of freedom is given by the m-theorem, proved 
in Ref. [95] : the theorem gives sufficient conditions for a loop integral to vanish in the large m -limit. For 
the one-loop case it concerns 

I = m a [d 4 q - fff (151) 

where 

N 

k = q + XI ^'j Pj' m < = or m ' ( 152 ) 

P(q) is a monomial in the components of q, {p} are the external momenta and a is an arbitrary real number. 
Let co be the IR degree of / at zero external momenta; we define 

d = diml, ^ = min{0,co}. (153) 

If / is both UV and IR convergent and d <Q. then / — > when m — > oo. 
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In conclusion one should say that BSM Lagrangians can be also classified according to decoupling. 
Thus the strategy can be summarized as follows: first, fix benchmark models to parametrize deviations from 
the SM, then search for 



benchmark models 



d = 6 operators G {J^dec} 

G {-^non-dec} 



6.6 Mixing 



There is one assumption in Eq.(15) and in its interpretation in terms of ultraviolet completions: the 
absence of mass mixing of the new heavy scalars with the SM Higgs doublet. Presence of mixings changes 
the scenario; consider for instance a model with two doublets and 7 = 1/2 (THDM), and §2- These 
doublets are first rotated, with an angle j8, to the Georgi-Higgs basis and successively a mixing-angle a 
diagonalizes the mass matrix for the CP-even states, h and H. The SM-like Higgs boson is denoted by h 
while the VEV of H is zero. The couplings of h to SM particles are almost the same of a SM Higgs boson 
with the same mass (at LO) only if we assume sin(j3 — a) = 1. Therefore, interpreting large deviations in 
the couplings within a THDM should be done only after relaxing this assumption. 

The case of triplet-like scalars is evem more complex; in the simplest case of a triplet with Y = 1 there 
are four mixing angles, all of them entering the coefficients of 



1 



ts-/(t s ) 



(154) 



in the amplitude for h — > yy (where S = H + ,H ++ ) and giving the couplings hH + H and hH ++ H , where 
h is the SM-like Higgs boson. Only in a very special case, requiring also zero VEV for the triplet, these 
couplings assume the simplified form 

M 2 M 2 

c hH+H- =2 -^-> c hH++H— = 2 — — > ( 155 ) 

where v is the SM Higgs VEV. Furthermore, decoupling of the charged Higgs partners depends on the 
mixing angles and it is the exception not the rule. 



7 Decays into 4 -fermions 



With a light Higgs boson the decay H — > VV is not open, and one should consider the full H — > 4f 
channel. In order to understand how the calculation can be organized we start with H — > ZZ where both Zs 
are real and on-shell. 

• H^ZZ 



The SM amplitude is 

_ _ M r r SM;L0 g 2 SN 



8^ V + 



16k 



2'T 



where 



v M 
Pi ' Pi 



(156) 



(157) 
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We introduce auxiliary coefficients 



A*=A 5 K s e ±Alc e , 



4 - 



Al+Al+IA 



5K- 



The full amplitude reads as follows 



86 M 2 
4V2 c e 
M 2 

c 2 



v = 2 5 / 4 Gp /2 (^ D + Jh , 



(158) 
(159) 



56 



— 4 

G^M ,^SM,NLO 



2V2c 2 n 2 D 



4y/2 
1 



8A^s e c e -A| 



^M 2 
V2 H 



g6 

4^2 



8Ay s e C0 

AyCQSQ — Ay^ — AyCg — — A R — 

4 ce 



-A° K ) 



(160) 



GfM ^-,sm,nlo 

2y/2&n 2 T 



1 



g6 

4^ 



8A v f e c e -A«) 



(161) 



Following the same strategy we consider the case M H < 2M Z . The process to consider is then 
• H^Zff 

which means H — > ZZ* — > Zff. If we work at LO and only include local operators proportional to the 
tree-level HZZ coupling, then the SM amplitude is multiplied by a factor 



r L z 



/SM 



1 



86 K Z 

'■ ^AC 



4^2 



where the correction w.r.t. the SM is given by 

- M 2 ^ 



k-:, = =^c e [c e Al-s e Ai) -A^ + AA dK -As e c e A\ 



Therefore we can use (at LO), the SM result and write 



r(H -> zz*) = £ r (h -> zz* -> zff) = (l + ^= K-f c ) r SM (h -> zz*) , 



/ 

where the SM partial width is 

r SM (H^ZZ* 



647T 3 H Im 2 



40 2 190^ 



(162) 



(163) 



(164) 



(165) 



F is the three-body decay phase-space integral, 



Fix) = (x - 1) (jx - H + I) + | (1 - 6x + 4x 2 ) In; 



„ 1 -8x + 20;t 2 
+ 3 ; — arccos 

Note that this result cannot be extended beyond LO. 



1 3x- 1 

2^72" 



(166) 



28 



• II wiT 

Similarly to the previous case we have a correction factor 

M 2 



M 

and the partial decay width can we written as follows 



K: c = =^Al-A^ + 2Ai + 4A dK , (167) 



3GlM 4 — ( M' 



r SM (H^ WW*) = ^3-M H F [ — . (168) 



327T 3 \Mh 



Taking the ratio we obtain 

K/.w — — rrr; ^zw ( 1 -r-^v ) (169) 
where the correction factor is 

r zw = 2A| + 4 ^c fl A^ + =| f e (i e A 4 K + c e A|) . (170) 

• H^ZZ^ffff 



r(H^zz*) / g6 
r(H^wr) zw lv 2V2' 



However, if we want to deal with the whole process without approximations the final state is 4-fermions 
(say ee|xu) and the SM amplitude has also non-factorizable contributions (e.g. pentagons). 

-#sm = ^fc V (Pl,P2) A Ma (pi) A VJ g (p 2 ) J P {qiM) +^nfc(Pl,P2) , (171) 

where / is the fermionic current 

J^(q,k)=gu(q)f (v f + a f f)v(k), p = q + k. (172) 

Furthermore, A^ v (p) is the Z propagator and ^# n f c collects all diagrams that are not doubly (Z) resonant. 
The question is: can we extract informations on 

2 2 
77SM 77SM,LO 1 S j-iSM.NLO r SM 5 j-^SMjNLO / 1T! i 

Fd =F + 16^ Fd Ft = 16^ Ft ' (173) 

or deviations from the two SM structures from the decay H — > 4f? 

The form of the Z propagator depends on the choice of gauge but, as long as the fermion current is 
conserved all differences are irrelevant. With the polarization vectors of Appendix B one obtains 

£ e± tl {p,X)e* Lv {p,X) = 8^ lv -^^--e Ltl (p)el v (p). (174) 
x=-i,+i P 

and we can safely replace the 8^ v in the propagator with a sum over polarizations, even for off-shell Zs. 
Using Eq.(174) we replace 

A"» -+X> M (p,A,K(p,A,)A(p 2 ), A(/) = — \=- 2 , (175) 
x s + M o 
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with M 2 , = M 2 /c 2 e , p 2 = -s and 



en(p,0) = e£(p), « M (p,±l) =«^(p,±l). 
We introduce the following matrices 



ep(pi,i)e v (p 2 ,j), 



Dij(p) = I Ei(p)E](p), Eiip) = J" (q,k) e* v (p,i) 



(176) 

(177) 
(178) 



spin 

where i, j = — 1 , 0, + 1 and p = q + k. We obtain 

2 



I 

spin 



I^4A*Cpi)£»,7CP2) A( Si )A(j 2 ) = J>y« A( Sl )A(s 2 ) 

ijkl ijkl 

-i 2 

A UU + L A W + L A tf«J A 0*l ) A (S2) 

ij Ktfi 

¥J 



(179) 



where ^# is the matrix element comprising all factorizable contributions, not only the SM ones. A,,,-; gives 
informations on H decaying into two Z of the same helicity (0,0 etc.), Aipj on mixed helicities (0, 1 etc.) 
while the third term gives the interference. Therefore 



A; — A ,7,7, 



A.. — A.... 



(180) 



are good candidates to define pseudo-observables. The final step is achieved through the realization that 
pseudo-observables are defined in one-point of phase-space and the choice must respect gauge invari- 
ance [96]. The amplitude in Eq.(179) has the general structure 



Y j aij{s,s\,s 2 , ...) A(ji)A(j 2 ) 

ij 

= Y, a ij ( s b,s z ,sz •••) A(si)A(s 2 )+N(s,si,s 2 , ...) , 

ij 

where N denotes the remainder of the double expansion around s\ i2 = sz, s = — (p\ + p 2 ) 2 and 

1 



(181) 



A(s) 



s-s z 

sh,sz being the H,Z complex poles. Therefore, we define pseudo-observables 



r,- = / d<£>i^ 4 a a {sn,sz,sz ■ ■ •) A(ji) A(s 2 ) 

spin 



(182) 



(183) 



with similar definitions for r,y. Since the problem is extracting pseudo-observables, analytic continuation is 
performed only after integration over all variables but s\ ,s 2 . Nevertheless, if one wants to introduce cuts on 
differential distributions alternative algorithms must be introduced, see Ref. [97]. 
The matrices D,E are given by: 



Poo = ~ (sis 2 ) 1/2 z H 



,-4 



SlS2 yy 
4 
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P++=P—=-i 
1 



1/2 



sis 2 N[Nl 

r2\"l/2 



d-^t) £(^1,^2,91,92), 



'd -^t) |2jiJ 2 *45 



^46^56 — *6 ZH 



ZH 



}■ 



where we have introduced 

Sy= Jj + Jj, Zh=^h-^i-^2, £(^1,^2,91,92) = %v«/3^f ^29? 9 2 - 
The elements of the D -matrix are given by 

D 00 (pi ) = 2 (V 2 + V 2) [4 + _L (2,! , 2 - (4 + 4)) 

D (pi ) = D ++ = (Vl + V 2 ) ^ { (* 34 + *5 6 ) *56 

+ \^[ 2(s 2 l sl-(sl 4 -2s l s 2 )s 2 56 )-(s 3 4-s 56 ) 2 s l s 2 -(s 3 4 + s 56 ) 2 s 2 56 }] 

D -(P!) = 7 ^ w-^) [2^4fe4-, 56 ) 2 ; 



- (^+ + V-) {^56 + ^ [(S34 + S56) (S! S 2 - 4) - 2 (4 - Sl S 2 )s 56 \ } 

7T ( V + ~ V -) (^34 - S56) (^56 ^34 ~ * 1 S 2 ) 



-S1S2) + 



1 1 



(534 - s 56 ) 2 si s 2 + (534 + s 56 ) 2 4 



(184) 
(185) 



(186) 



2Ayvj 

+ W + v?)^{-(4.....„. 4ML 

- 2 (s 2 s 2 . - (4 - 2*i J 2 ) 4)] } 
and similarly for Dij(p 2 ). 

The result of Eq.(179) does not include non-factorizable diagrams. To include them we will follow the 
work of Ref. [98] where standard matrix elements (SME) are introduced (see Eq. (3.1) and Eq. (3.2) of 
Ref. [98]); they are made of products of 

^° = lu(qi)rn {l+o^)v(ki), r i ^ va = ^u(q i )Y^ v Y a (l + oy 5 ^-), (187) 
with = ±1 and i = 1,2. For example one has 

^12,« = r l,«; Mr £x a? V£« (18g) 

etc. The non-factorizable amplitude becomes a sum 



(189) 



where the F are Lorentz invariant form-factors computed up to NLO but excluding those that are double 
resonant; the full answer follows by adding this amplitude to J%f c . Note that 

(q h ki) =g (vf + Of) T^+g (vf-Of) rji". (190) 
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8 Double Higgs production 



A non-zero value of a g fives a contribution also to the ggHH vertex, contributing to double Higgs 
production, gg — > HH (see also Ref. [99]). 



]Xa 

X 


ggHH 


4a g G Fg6 T flv 5 a ' h . 


(191) 


vb 









An additional contribution to double-Higgs production come from the HHH vertex where py,p2,P3 are 
the momenta of the outgoing bosons with p\ +p2 +P3 = 0. There are also quartic couplings 





9 Perturbative unitarity 

In this section we study constraints from perturbative unitarity. With no informations on the Higgs boson 
mass there are two different scenarios in V L V L — > V L V L scattering: 

1. M^,Ml<^Ml<^s 

2. M^,Ml^s^Mi 

Assuming a light Higgs boson we analyze a new option, 
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The SM result iw well-known, given 

d 



^v L v L ^v L v L 



T(s,t) 



we derive 



with a critical mass 



C(w+w L ^w+w L 



L0 ~ 16 



1_ yo 



A r, , 



5 — ). oo 



r° (m h =m c 



i. 



(194) 



(195) 



(196) 



Anomalous couplings violates perturbative unitarity. However, one has to be careful in formulating the 
problem: the region of interest is 

• M^M^M 2 <s<A 2 . 

When s approaches A 2 the effective theory must be replaced by the complete renormalizable, unitary La- 
grangian and it makes no sense to study the limit s — > oo in the effective theory (for a discussion see 
Ref. [100]). To summarize, anomalous vertices with ad hoc (scale-dependent) form-factors are frequently 
used but one should remember that they cannot be put down to an effective Lagrangian. 

However, it is well known that heavy degrees of freedom may induce effects of delayed unitatity can- 
cellation in the intermediate region and these effects could easily be detectable [101]. Without using the 
equivalence theorem, we compute 



[ SM+AC 



1 



167TX(s,M 2 ,M 2 ) J- S+ 4M 2 



-t s 



dt r SM+AC (w+W L -> W+W L ) , (197) 



with a cut to » M 2 /s to avoid the Coulomb pole. Longitudinal polarization vectors are defined as fol- 
lows [102,103] 



e]i{P\) 



2 2 

{pi -P2Pln+M 2 p2n) e^(p 2 ) = — — (pi -P2P2h+ M2 PIh) 



Ms/3; 



M 



2 2 

e \i(Pl) = T7 o (P3-P4PlH+ U2 P4n) «J;(P4) = =-^~ {P3-P4P4n+M 2 p 3lJ ,) , 

MjpM M^Pm 

with = 1 - 4M 2 /s. In the limit Mi^,M\,M\ <Cs<A 2 we obtain the following result 



(198) 



1 SM+AC 



l - {2 + 5t -ti) (l-to)c e $e (l - 2l|) 2^4; 



+ (!- ? c)) 2 U K + Ai-A dK -6 S f s A 



+ [i (1 1 + \0t - 13 tl) ) - 2 (1 + 2 to - 2 tl) §1 ] c e s e 4 } ^ # 6 



+ 



16^2 



(l-'o) 2 ^A|-A 4 K 



7T 



(199) 



As expected the SM part contributes to the constant part while the part proportional to has positive powers 
of s (up to power two). The leading behavior is controlled by the &\ operator. 
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10 Conclusions 



We have described possible deviations from the Standard Model parametrized in terms of effective d = 6 
operators made of Higgs, gauge and fermion fields, without making the hypothesis that the new physics 
shows up in the Higgs sector. Furthermore, we allow effective operators generated at tree level and by loops 
of heavy particles. 

In this paper we have discussed the implementation of effective Lagrangians with emphasis on renormal- 
ization. Examples of Lagrangians producing the d = 6 operators have been shown and we have discussed 
both the decoupling and the non-decoupling scenarios. In agreement with the work of Ref. [100] we have 
been following the effective field theory approach which is cleaner than that of anomalous couplings. An 
effective field theory is the low-energy approximation (E <C A) to the new physics and it is only useful up 
to E « A: above A it should be replaced by a new effective theory, parametrizing the low-energy effects at 
a yet higher scale. 

Effective theories should not be considered beyond their UV cutoff, although this is often done in the 
literature with the introduction of methods for unitarizing the model, e.g. form-factors are introduced; this 
requires specific assumptions and cannot be formulated in terms of an effective Lagrangian. 

There are many scenarios, e.g. an interesting one (see Ref. [87]) has no new charged fermions and only 
new bosons. This would unambiguously rule out a large class of BSM theories. There are also scenarios with 
new physics which will be extremely difficult to distinguish from minimal SM, e.g. see Ref. [104]. However, 
the analysis of all possible options should not be done hiding uncertainties or the bias from discovering using 
the minimum p -value. Opportunities for precision measurements and BSM sensitivity have been recently 
described in Ref. [105] and in Ref. [106]. 

One final comment is needed: the strategy described in the Introduction amounts to search for deviations 
around a minimum which we assume to be SM. If measured deviations will be large, we will face a problem 
of interpretation: indeed, consider the ratio 

R = ^L, (200) 

where ghvv is the tree-level coupling of the scalar resonance h to VV; if h = H, the SM Higgs boson, then 
R = l. Assume that R sxp turns out to be close to —1/2, this will be hard to interpret in terms of a weakly 
coupled theory and it becomes questionable to trust predictions from an effective Lagrangian, based on 
d = 6 operators, with Wilson coefficients of that size; to state it differently, d = 8 operators and insertion of 
d = 6 operators in SM loops are all equally important. However, some anomalous value of R exp could very 
well be close to another weakly coupled theory; for instance, the h 5 of the Georgi - Machaceck model [80] 
has R = — 1/2. Starting from the new weakly-coupled Lagrangian will allow us to trust the prediction. Of 
course, one would like to be as model-independent as possible without repeating the fit for many different 
starting points; however, there are only very few representations of SU (2)l <8> SU (2)r that respect custodial 
symmetry, and they should be included in a more comprehensive analysis. 

A recent note by ATLAS Collaboration [107], using data taken in 2011 and 2012, reports that, within the 
current statistical uncertainties, no significant deviations from the Standard Model couplings are observed. 
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A Appendix: The ghost Lagrangian 

In this Appendix we give the explicit expression for the Faddeev-Popov ghost Lagrangian. 

_ _ _ M 2 — 

^f FP = X-d 2 X- +X + d 2 X + + Y z d 2 Y z + Y A d 2 Y A -M 2 (X+X+ + X-X") Y Z Y Z 

c e 

_^m[+(4 + 2^ k ) ^t] (X + X+ + X-X-)H 



+ l - 8 M (i+4^r) (x+x + -x-x-)^< 



1 M 

2 c 2 



l - g -^ + (Aa\c e se-a^) (y z X"4> + - Y z X+<)-) 



+ 



i M 



1 + [SdySece +a K + 2a dK 
M 2 



M 2 
A 2 " 



Y Z Y Z H 



~2 ~2 

^e ~ c e 
ce 



+ 



Aa\se (1+24) +4 



3 (4 + 3c 2 ,iM 2 



-^}(X + Y Z 4> + -X-Y Z *-) 



+ igs e [\ + (Aa\s e +alce) ^jr] (x~ d^X~ -X + d^X + )A^ - (X+X+-X-X") 



1 + [Aa\c e s e +a| 



\ M 2 


ce" 




/ A 2 " 


*e J 


( 


\ M 2 


r 2 i 

c e 




/ A 2 " 


V 2 





Y A d,X+-X-d,Y A ) W" 



+ (y a x + -x- y a ) a M w-- (y a x--x+y a ) d M w+- (y a ^x--x+^y a ) w+ 



+ igMs e 
+ igce 



1 + (4ay<?ese +a K 
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M 2 c 2 
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+ ^ce 



M 2 
A 2 " 
M 2 



(x - d M x- - X + d M X+) z M - (X +X+ - X - X") d M 



l-(44c e f e +4) [(Y Z ^X + -X-^Y Z ) W m + (y z X + -X-Y z )^W^ 



(y z x--x+y z )^w+-(y z ^ lVzr . : 



,X--X + d u Y 7 W 
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B Appendix: Polarization vectors 

A convenient choice for the polarizations in H — >■ VV is the following: 

ehn(pi) = -N[(pi ■ P2Pifi+sip 2 ^i), e Lfl (p 2 ) = -Nl (pi ■ p 2 p 2il + s 2 pin), 
where N\ i2 are the normalizations, pf = —57 and 

e±n{pi,^> = [np(pi) + iXNn(pi)], N^(pi) = (siy l/2 e^ a p p n a (p^e^pi) p p , 

in (pi) = iN[ e^pp kf p\ p p 2 , n M (p 2 ) = iN\ k% p\ p p . 
With this choice one obtains 

£ «±M(^'^) g iv(P>^) = 5 /iv-^T^-«L/xO , ) e Lv(p)- 
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Using P — > p\ + p2 2 — > q\ + k\ + q 2 + k 2 we define 



P\-P\=S\ P2-P2 = S2 Pl-P2 = ^{s]_+S 2 -S] i ) 

q\-qi=0 q 2 -q 2 =0 k\-k\=§ k 2 -k 2 =0 
1 



qi-ki = 
k\-q 2 



-s\ q\-q 2 



1 , 1 

-53 q\-k 2 = --SA 



~s 5 h-k 2 = -^s 6 q 2 -k 2 = -^s 2 , 



where we allow for an off-shell Higgs boson, P 2 = —su- We derive 



-1/2 



1 



n l = 2^(su,si,s 2 ) , 



where X is the Kallen function. Furthermore, 

1 



N 



i_l 



N, 



1 



2_L 



(s 5 +s 6 ) (s 3 + S 4 ) ~SiS 2 
(s 4 +s 6 ) (s 3 + S 5 ) ~SiS 2 
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